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ABSTRACT 



The effects of cutouts on the strength of glass 
reinforced plastic <GRP) for use in naval ship hulls have 
been studied. Experiments were conducted on an E-glass 
fiber/polyester resin laminate provided by the U.S. Navy 
Minesweeper Hunter <MSH> shipbuilder in 48 in. x 15 in. 
panels. Tensile tests were run on full-sized panels with 
circular cutouts of 3.75 to 7.50 inch diameter with and 
without edge notches; conventional tensile and fracture 
toughness tests were also run. 

Three theoretical models were evaluated. First, the 
notch sensitivity was evaluated, and the material was found 
to be notch sensitve, so that models which include stress 
concentration effects are required. A simple stress 
concentration model appears valid, but can only provide an 
accurate prediction when the measured tensile strength is 
increased by a factor to represent the actual local strength 
at the cutout edge. Linear elastic fracture mechanics 
provides strength predictions reasonably close to the 
experimental results when cracks are present at the edge of 
the cutout. 

In addition to the direct experimental results of the 
study, the effects of a seawater environment and cyclic 
fatigue loading are estimated based on earlier studies. The 
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overall effects on ship hull design safety factors are 
calculated. The required safety factors are excessively 
high unless cutouts are reinforced to eliminate the stress 
concentration effect. Detailed recommendations for further 
research and development needed to advance failure analysis 
in GRP hulls are provided.' 
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Title: Professor of Polymer Engineering and 
Professor of Civil Engineering, and 
Principal Research Associate 
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CHAPTER 1 



INTRODUCTION 



This thesis investigates the effect of flaws on Glass 
Reinforced Plastic (GRP) strength. Experiments have been 
designed to measure the effect which holes with edge cracks 
have on laminate strength in a GRP plate of E-glass fiber 
(woven roving and chopped strand mat) in a polyester resin, 
which is a typical naval ship hull material. This material 
is of particular interest because of the impending 
construction of the U. S. Navy's 434 ton Minesweeper Hunter 
(MSH) , which is under contract to Bell Aerospace 
Corporation. This represents the U. S. Navy's first attempt 
to use GRP in a large ship hull, thus the program is under 
scrutiny both from within and outside of the Department of 
the Navy. 

This thesis presents experimental data which, when 
correlated with the appropriate theoretical failure models, 
may provide a better understanding of the relationship 
between flaws and strength in a GRP naval ship hull. It is 
hoped that any new understanding will permit more reasonable 
factors of safety in the design of GRP ships which should 
allow GRP to be considered for the new generation of 
weight-limited, high performance ships required by a modern 
U. S. Navy 123 . 
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CHAPTER 2 



BACKGROUND 



Qverv i ew 

The use of glass rein-forced plastics in the ocean 
engineering industry began in World War II. At that time 
GRP represented a new material with the strength, stability, 
low weight, environmental resistance and ease o-f fabrication 
needed for high-volume production of large marine parts such 
as boat hulls C33. Only the high cost of raw materials and 
the slow, expensive processing methods of this period 
prevented its widespread use as a hull material. As shown 
in Table 1, the advantages of GRP now outweigh the 
disadvantages in naval ship hull and structural uses C4] . 
The significance of the advantages are clear when we view 
the demand for GRP during the last ten years. Table 2 
presents data for GRP end use and shows that the growing 
total demand has been paralleled in the marine industry C5D. 
At the present time GRP materials are being used in the 
construction of boats, buoys, control surfaces, deckhouses, 
fairings, masts, propeller shrouds, protective coatings, 
radomes, sonar domes, pipes and other applications. 
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Hul 1 Struc tures 



Of particular interest is the use of GRP in naval ship 
hulls. The U. S. Navy pioneered fiberglass boat 
construction with the production of a 28 foot hull in 1947 
and was influential in the use of hand lay-up procedures 
using woven roving and fire retardant polyester resins for 
boat hulls [63. In 1951 a development contract was awarded 
for the construction of a minesweeping boat (MSB) with a 
structural sandwich core hull. The 57 foot boat was 
launched in 1956 and served satisfactorily for over 10 years 
[7,83. By 1973 the Royal Navy had launched what was then 
the world's largest fiberglass ship, H.M.S. WILTON, a 153 
foot, 450 ton minesweeper which used a mono-skin for its 
hull [93. The success achieved in the H.M.S. WILTON 
prompted most of the world's navies to adopt the use of GRP 
in mine warfare ship hulls. Table 3 is a summary of 
existing mine warfare ships [103. However, these ships have 
not been without their problems, the most predominant may be 
the presence of defects and their effect on the structural 
integrity of the ship hull [113. 
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Hul 1 Defects 



Thomas and Cable [11] have investigated the important 
defects in GRP naval ship hulls and classify the driving 
material defects as Type A or Type B. Type A are gross 
defects while Type B are more subtle and include voids, 
de 1 am i nat i ons , reinforcement discontinuities, holes and 

notches. The Type B defects are considered to be more 
common and to have the greatest potential for property 
degradation [12]. 

With the use of modern Quality Assurance procedures all 
of the Type A and many of the Type B defects can be 

detected, evaluated, and corrected during production. 
However, when the remaining defects are placed in the 
context of a naval ship hull, the most deleterious effect 
was predicted to be [11] due to unreinforced hull cutouts 
(holes) and the resulting stress concentrations. This 

conclusion is based upon the tremendous number of relatively 
large hull penetrations that must exist in a ship for 
access, piping and wiring runs. The magnitude of the 

problem increases because these openings will have to be cut 
out of the finished laminate and may range in size from a 
1/4 inch diameter wiring run cut to a 5 foot x 10 foot or 
larger opening in the main deck of the vessel. The larger 
cutouts would presumably be reinforced. However, the worst 
problem may arise because hull penetrations cannot be 
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fabricated totally crack free. The effect that small edge 
cracks have upon the engineering properties of the GRP 
laminate and the overall ship design elements may be a major 
area of concern, particularly considering the relatively low 
fracture toughness of typical GRP ship hull material as 
compared to ship steels. 

The remaining chapters of this thesis are dedicated to 
understanding the various theoretical failure models and the 
correlation between experimental data and these theories. 
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CHAPTER 3 



THEORETICAL MODELS 
OF FAILURE FOR GRP 

In order to estimate the effect that a hull penetration 
and its associated cracks has on GRP laminate strength, 
several models which predict the performance degradation 
were evaluated. Each o-f the models addresses a 
through- the-th i ckness penetration in a 6 ply GRP laminate 
described in Table 4. The assumptions used in each model 
along with the limitations are presented. Detailed model 
development and calculations are provided in the 
corresponding appendices. 

Notch Insensitivity Model 

This model deals with a quas i -due t i 1 e material where a 
notch or crack is of no consequence in failure, beyond the 
reduction in cross-sectional area. In composites this 
behavior may result from widespread matrix cracking (damage 
zone) and crack tip blunting that occur at a crack site 
1 1 3 3 . The remaining strength of the specimen, is 

approximately the load averaged over the remaining cross 
sectional area. For the geometry presented in Figure 1, 
failure would occur at 

^ net = E 

t <w - 2a - D) 
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where P is the load at failure. In terms of the material 



tensile strength, u ts> the behavior is notch insensitive 
if 



ne t “ 




< 2 ) 



Appendix A presents calculations for various hole and crack 
sizes in a panel of typical GRP naval ship hull material 
using this model as the basis for calculations. 

Stress Concentration Model 

In this case failure is predicted by using the elastic 
stress concentration factor, K T , for the hole, without 
considering edge cracks where present. Kj is defined as 
the ratio of the maximum stress in the presence of a 
geometric irregularity or discontinuity to the stress which 
would exist at the point if the irregularity were not 
present [143. This relates to the failure stress as 

<=r_ failure ~ ^^^"uts ^ ^Tor thotrop i c 

In order to apply this model, the stress concentration 
present in an isotropic, finite width plate presented in 
Figure 2 is first calculated. Savin [153 has tabulated 
values of <a e/P along the contour of a circular hole for 
various values of X= a/b. As expected, the most dangerous 
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stress concentration exists at the edge where 6 



90o , 



Appendix B presents calculations for the maximum isotropic 
stress concentration factors for various hole sizes as well 
as the expected or predicted failure stress for a typical 
GRP naval ship hull material. 

The typical naval ship hull material is orthotropic, 
and so the effect of an i sotrop i c i ty cannot be ignored. A 
correction factor, ^ , can be calculated by relating the 
orthotropic stress concentration factor for an infinite 
plate to the isotropic stress concentration factor for the 
same infinite plate. For the orthotropic infinite plate in 
Figure 3, the stress concentration factor is given by [163: 



Since the circular hole is a special case where R=a=b, the 
stress concentration factor for an infinite plate with a 
circular hole becomes 



*^Tor thotrop i c — 1 + n <a/b) 



< 4 ) 



^Tor thotrop i c — 1 + n 



(5) 



where 



n 
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and E x f E y , G xy , and V^ x are the orthotropic elastic 

constants. For an isotropic material, n=2 and equation 

< 4) becomes 1 1 63 : 

^Tisotropic =1+2 <a/b> <?) 

and -for the case with a circular hole 

^Tisotropic = 3 <8> 

The effect that anisotropy has on the stress concentration 
•for the geometry presented in Figure 3 can now be 

calculated. The correction factor is used to correct the 

stress concentration factors presented in Appendix B. 
Appendix C contains the calculations for the correction 
factor and Table 5 presents the corrected stress 

concentration factor and predicted failure stress for the 
orthotropic material. The factor for the material in 

question gives a stress concentration of about 15% higher 
than the isotropic case. 

Linear Elastic Fracture Mechanics Model 

In this model the material fails in a brittle fashion 
and the classical linear elastic fracture mechanics approach 
prevails. Fracture mechanics can be used to predict either 
the stress < *3"" ) or crack size <a> required to fracture a 
particular sample. Broek [133 relates these two parameters 
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•for the -finite plate in Figure 4, with the stress- i n tens i ty 
■factor Kj suc h that 



where w is the plate width. 

Kj is a -function of both specimen and crack geometry. 
In the case under consideration, the hull panel will be 
modeled as a flat rectangular sheet containing two edge 

cracks at a circular hole which is then subjected to uniform 
uniaxial tensile stress. This geometry is depicted in 
Figure 1. Given this geometry and loading scheme, a 

compendium of stress- i n tens i ty factors can be consulted for 
the Kj relationship [173. Figure 5 presents the Kj/K q 
relationship versus a/b for various R/b ratios and is 
accurate to within V/.. The Kj/K q relation is used to 

correct the K Q value for the geometry considered by using it 

i n 




<?) 



*1 = K 0 <K,/K 0 > 



< 10 > 



where 



K q = <3 



< 1 1 ) 
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Formula (10) is -for an isotropic material and the subject 
material is an orthotropic GRP. However, earlier work has 
shown that orthotropic values differ by less than 25X [183 
and the isotropic value is used here in the absence of an 
orthotropic analysis o-f the problem. 

Failure, in the LEFM case, is defined to occur when the 
predicted stress- i n tens i ty -factor <Kj> greater than or 

equal to the critical stress-intensity -factor Kq , the 
material fracture toughness. The subscript Q used with Kq 
indicates the candidate opening mode critical stress 
intensity factor as suggested by ASTIi (E399), and is 
distinct from Kj c> as the latter satisfies certain validity 
requirements which are not available yet for composites 
[193. A Kq value of 17.5 ks i < i n > j s expected based upon 
K-cal i brat i on and test data for a similar E-glass fiber and 
polyester resin material [183. However, a new value for the 
particular material in question has been determined using 
the test geometry in Reference 18. 
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CHAPTER 4 



MATERIAL AND 
TEST METHODS 



Mater i al 

In order to evaluate each of the failure models and 
select the one that best represents actual material failure, 
numerous tests were conducted on a candidate material. The 
material selection process was simplified by the fact that 
the only GRP naval ship hull under construction in the 
United States was the MSH. The test samples selected were 
actual panels of this hull material. The test panels were 
fabricated at Bell Aerospace Corporation's MSH production 
facility in New Orleans, Louisiana and shipped to MIT for 
analysis. The panel dimensions were 48 inches x 15 inches x 
6 ply (.3125 in. nominal thickness) and were designed to 
have the material properties reported in Table 4. 

Test Methods 

The first test method concerned the effect of holes, 
both with and without notches, on the strength of large GRP 
panels. A series of 16 panels, 48 in. x 15 in., each with a 
central cutout, were tested under tensile load to failure. 
Ten of the panels had holes without edge cracks and 6 of the 
panels had holes with machined edge cracks. Table 6 
tabulates the test specimen characteristics and the method 
used to cut the central hole. 
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Specimens were loaded in the specially designed and 
•fabricated grip system shown in Figure 6. The grip was 
designed to trans-fer the tensile load to the specimen using 
pin loading. Each grip consisted o-f two 3/4- inch thick 4340 
heat treated steel plates. Each plate was bored to receive 
a 1-1/8 inch pin and -f if teen-5/8 inch grade 8 bolts through 
which the load was transferred to the specimen. 

The GRP samples were drilled with a carbide tipped 
drill to enable them to accept the fifteen-5/8 inch bolts. 
An aluminum drill jig with RC 60 drill bushings provided the 
means by which precise hole alignment was maintained. Each 
specimen was then fit between the plates, bolted together 
and torqued to 50 ft-lbf. Figure 6 presents the plan for 
the grip plates and Figures 7 and 8 show a specimen bolted 
in the test apparatus. Figure 9 is the 600,000 lb electro 
hydraulic screw machine at the U. S. Army Materials 
Technology Laboratory in Watertown, Massachusetts that was 
used in testing the samples. 

A chief concern in this test was obtaining a uniform 
load distribution across the entire sample width. By 
design, the applied load at the 1-1/8 inch pin was required 
to be transferred to the sample by the 15 bolts. This was 
verified by monitoring the strains at four equidistant 
points across the width of the sample. Strains up to 0.5/ 
were measured and an even strain distribution was observed 
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across the sample. Figure 10 is a load-strain curve For the 
four strain gages. 

The second test method utilized the same grip system 
but measured the ultimate tensile strength o-f a large 
specimen. Three -full size test panels were machined to the 
specimen shape indicated in Figure 11 using a fine-toothed 
diamond bladed band saw. The shape of this specimen was 
determined to give the highest stress at failure without 
breaking in the grips C193. In this test as well as the 
first a <600,000 1 bf Baldwin e 1 ec trohydrau 1 i c screw type 

universal testing machine was used, with a load rate of 6000 
Ibf/minute. Laboratory temperature and humidity were used 
in all tests unless otherwise noted. 

The third test was a conventional tensile test using 
the specimen shape in Figure 12 C19J. Six samples, three 
cut vertically and three cut horizontally from full size 
panels, were machined to this shape using a diamond-edged 
wheel and a TensilKut router. Sample strain was measured 
during this test using an extensometer which allowed the 
bulk modulus of the samples to be determined. An Instron 
universal testing machine was used for this test, with a 
displacement rate of 0.05 inches/minute. 

Finally, a conventional fracture mechanics test method 
was employed to determine the fracture toughness of the 
material. Six specimens, three cut vertically and three 
horizontally from full size panels, were machined to the 
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shape in Figure 13, -from Reference 18. For the cleavage 
specimen the isotropic value of K Q is g j ven j n [20] as: 



Kn = 2^3 P <C + 0.6? > <12> 

Q iW h 

where 

P = applied load at fracture 
B = thickness of specimen 
H = half width of specimen 
C = crack length 

For the specimen shown in Figure 13, Equation (12) reduces 
to 

Kq = P ( 6 . 03C + 6.24) (13) 



The geometry used in this specimen permits this test to only 
approximate the stress intensity factor due to free end 
effects. Data from this type of test has been in good 
agreement with finite element data given by Kanninen E213 
for the specimen size and crack lengths of interest in this 
test. However, a more accurate value for Kq can be obtained 
using the K-c a 1 i br a t i on data for a very similar 
E-glass/polyester material reported in Reference 18. Table 
9 presents both values for Kq. 

The procedure for measuring the fracture toughness was 
to load the specimen until fracture occurred at the precut 
notch. The maximum load at fracture was recorded and 
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provided linearity was maintained in the load-de-flection 
curve until the -fracture point, this maximum load value can 
be used to calculate Kq . This procedure is repeated a 
number o-f times by increasing the de-flection in order to 
keep the crack moving. This provides a number o-f data 
points although the load to cause initial crack extension 
was used in calculating Kq -for each specimen. An Instron 
universal testing machine with a displacement rate o-f 0.05 
inches/minute was used in this test. 

Finally, the e-f-fect on strength due to hole machining 
method was evaluated. Two o-f the 5.0 inch specimen holes 
were machined with an ordinary arbor mounted hole saw while 
■four others were machined with a carbide tipped parting tool 
mounted in a -flycutter. 



27 



CHAPTER 5 



DISCUSSION OF RESULTS 

The results of the -full scale tensile tests and the 
characteristics of the 1? large panels are presented in 
Table 7. Tables 8 and 9 present the results of the 
conventional tensile and fracture toughness tests. These 
experimental data are used to evaluate each of the three 
theoretical failure models discussed in Chapter Three. From 
this evaluation appropriate models are selected as best 
representing actual material failure. 

Notch Insensitivity 

The notch insensitivity model can be evaluated by 
calculating the degree of notch sensitivity displayed by the 
GRP. A good measure of this is the ratio of the strength of 
the remaining material, < ^ _ net» from Equation <1), to the 
ultimate tensile strength of the material, <T" U ^ S . A value 
of 1.0 would indicate notch insensitivity, while decreasing 
values less than 1.0 would indicate an increasing degree of 
notch sensitivity. Figure 14 shows one of the failed panels 
complete with machined crack and 5.0 inch hole. The average 
values of <3 ne ^ -for each hole diameter, notch length and 
machining method tested were compared to the average tensile 
stress of the three full size tensile coupons. Figure 15 
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gives these data as a -function o-f hole diameter plus notch 
length and shows that the material is strongly notch 
sensitive in all cases. UJhen edge notches are present, the 
behavior is more notch sensitive. The data do not cover a 
sufficient range of specimen sizes to show the expected 
increase in notch sensitivity with hole or notch size. 
These data also show that the failure loads measured for the 
two machining methods varied by 5.67%, with the hole saw 
method exhibiting the higher load values. The difference is 
not great enough to draw a definitive conclusion. However, 
it suggests that the machining method does not have a great 
effect on strength. These results indicate that the notch 
sensitivity model is inappropriate as a model for failure in 
this material, and would be very non-conservative. 

Stress Concentration 

The applicability of the stress concentration model can 
be evaluated by comparing the actual failure stresses 
measured in the test specimens to the predicted failure 
stress from the model. The actual failure stress in the 
specimens is taken as the far field stress at the failure 
1 oad as 

failure = P <14> 

t x W 

where P = failure load 



2 ? 



t = nominal specimen thickness o f .3215 in. 

Ul = nominal specimen width of 15 in. 

The predicted failure stresses for the hole diameters used 
in the tests have already been reported in Table 5. Figure 
16 shows a failed panel with a 5.0 inch diameter hole. 

A reasonable measure of the validity of the stress 
concentration model is the ratio of the predicted failure 
stress to the actual failure stress, with the model becoming 
more accurate as the ratio approaches 1.0. Table 10 
presents the data by which this model may be evaluated. 
When the raw data were compared, this model seemed to be 
i nappropr i ate for predicting the failure strength due to the 
low values of the ratio. However, a more careful review of 

the data was pursued. Figure 17 gives the ratio of 
predicted to actual failure stress against hole diameter and 
shows that for holes greater than 10 inches in diameter the 
model approaches validity, following a consistent trend. 
This trend may be explained by the fact that in very large 
holes the volume of highly stressed material more closely 
approaches the volume of stressed material present in 
ultimate tensile strength specimens. This is significant 
because the value from the ultimate tensile strength 
specimen is used in the stress concentration prediction, 
therefore, as the stressed material volumes converge the 
model becomes more accurate. 
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Figure 18 shows both the predicted and exper imen tal 1 y 
measured -failure stresses plotted against hole diameter. 
The plot shows that the predicted and experimental trends 
are similar, but that the predicted strength based on the 
ultimate tensile strength tests is lower. The difference in 
strength values may be explained by a magnification factor, 
which implies that the in-si tu ultimate tensile strength at 
the hole is higher than in a standard tensile test; an 
effect that is apparently due to the smaller material volume 
and the constraint that exists around the failure site, as 
compared with an ultimate tensile strength specimen with a 
narrow gage section (although the local strength at the hole 
cannot readily be measured experimentally). If this is the 
case, an increased ultimate tensile strength value is 
necessary before the predicted failure stress is in 
agreement with the experimental data. For the values 
reported in Table 10, the local strength must be 
approximately 1.58 times the measured ultimate tensile 
strength in the large tensile specimens. For the material 
tested, this would require an in-situ strength of 41,000 psi 
at the hole edge, a value which can be considered 
conceptually sound and of reasonable magnitude. As shown in 
Figure 18, the predicted trend with hole size using the 
higher ultimate tensile strength value is in good agreement 
with the experimental data for all hole sizes. 
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This model is not suf f i c i en 1 1 y accurate -for predicting 
the -failure strength o-f the material and geometry under 

consideration, based on the measured ultimate tensile 
strength. However, the stress concentration model does 

appear to be theoretically correct based on the trend o-f the 
data, i -f a more representative in-situ strength value is 
used. I-f the measured ultimate tensile strength is used, it 
represents a conservative method for predicting failure 
strength and may be of use in predicting failure for 
relatively large hole diameters. However, the method does 
not apply to either holes with edge cracks or small holes 

where the hole size is similar to the damage zone size, 

which is on the order of one inch for this material (Figures 
14, 16 and 19). 

Linear Elastic Fracture Mechanics 



The linear elastic fracture mechanics model can be 



eval uated 


by 


compar i ng 


the val ue 


of 


the stress 


i n tens i ty 


factor Kj 


at 


f a i lure 


predicted 


by 


the model 


wi th the 


exper i mental 1 y 


measured 


value of 


Kq. 


Using the 


techn i que 



described in Chapter three, the value of Kj for each unique 
crack length is calculated at the failure load. Appendix D 
presents a sample calculation and the Kj values for the 6 
panels with edge cracks that were tested. These values can 
now be compared with the experimentally derived value of Kq 
(using an orthotropic K-cal i brat i on) that has been reported 
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i n Tabl e 


9. 


Figure 19 shows a 


double-canti lever 


beam 


spec i men , 


comp 1 e te 


with crack 


propagation, used 


i n 


evaluating 


k q. 


The 


1 ongi tudi nal Kq 


values in Table 9 


are 


slightly higher 


than 


the calculated 


K j at -fail ure -for 


the 



holes with notches (Appendix D) . 

Figure 20 gives the results o-f the six test specimens 
reported in Table 9 and indicates that the value o-f load 
necessary -for crack extension is consistent with a constant 
Kq criterion using the analytic K-cal i brat i on o-f Re-ference 



18, over a 


broad 


range 


o-f 


crack 


lengths. 


The 


same 


K-cal i brat i on 


was 


used 


■for 


tests 


conducted 


i n 


both 



directions. However, they are only strictly valid -for 
elastic constants with a load in the longitudinal direction 
( parallel to the long side o-f the panel). Uith the 
experimental values -for Kq established, the applicability 
o-f the linear elastic -fracture mechanics model to holes with 
notches may be evaluated. The average Kq values were 25.9 
ksi(in)^' / 2 loaded longitudinally and 28.0 loaded 

transverse 1 y . 

Figure 21 compares the -failure stress predicted using 
Kq = 25.9 ksi ( in) 1/ 2 with the experimental data -for the 5 
inch diameter hole with di-f-ferent notch sizes. The 
predicted -failure stress is approached as the notch length 
becomes longer. This type o-f behavior is to be expected in 
this material because the damage zone size is o-f similar 
magnitude to the notch length, and precise modeling o-f short 
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cracks requires a two-parameter model C193. It should also 
be noted that the hole with notches calculation is based on 
isotropic properties, which could also lead to the small 
d i f f erences shown in Figure 21 C183. The significance of 

the model relies upon the fact that the predictions are 
reasonably close but non-conservative, falling above the 
experimental data. This model is considered valid but 
slightly non-conservative and therefore appropriate for 
predicting the failure strength of the GRP panels under 
investigation when edge notches are present or assumed at 
the cutout boundary, if the assumed Kq j s reduced slightly 
to provide a conservative prediction. 

Ship Hull Desion Strength 

The application of any model that predicts failure must 
be considered in conjunction with the other major factors 
which influence strength. In the case of GRP for naval ship 
hulls the predictions made using stress concentration or 
linear elastic fracture mechanics must be considered in the 
context of effects of both ship hull strength and fatigue in 
seawater . 

Fatigue studies for ships have chiefly concerned 
themselves with the maximum number of cycles a vessel will 
see at a particular stress level 1233. The area of most 
concern has been the maximum stress that a ship encounters 
in storm seas. Because of the lack of specific knowledge 
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concerning the size o-f storm waves to be expected at sea it 
has long been usual practice to design ships by assuming 
them to be poised on a stationary standard trochoidal wave 
C243. The standard wave height used by both Lloyds Register 
and the U. S. Navy is l.l(L)l/2 ( where L is the overall ship 
length. From this wave height a maximum bending moment and 
corresponding maximum ship structural stress can be 
calculated for a given hull geometry and thickness. This 
value is then subjected to a -factor o-f sa-fety -for the 
speci-fic material and load application used in the hull, and 
a maximum design stress -for the material is obtained. 
Thomas and Cable evaluated a ship hull structure similar to 
the MSH and report a maximum stress o-f 2970 psi, which 
corresponds to the maximum bending moment load on the 
standard wave [113. This conventional bending moment 

calculation ignores the dynamic e-f-fects o-f the ship's 
motion. The most deleterious motion results -from slamming, 
(the impact o-f the hull on the water during a large downward 
pitch). In naval vessels slamming stresses o-f a magnitude 
1.6 times the primary stress have been recorded C243. It is 
therefore reasonable to assume that the MSH type hull would 
see stresses 1.6 times 2970 psi or 4700 psi. 

Fatioue E-f-fects in Seawater 

The -fatigue o-f GRP in a seawater environment is a major 
concern. Investigators using a similar material to that 
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used here have found that the stress which can be withstood 



decreases by about 11. 6 /. of the ultimate tensile strength 
per decade of cycles when subjected to a sinusoidal 
load-time curve between a maximum tensile load and a minimum 
load which was approximately 5% of the maximum C26] . Figure 
22 reproduces these data. These results are applied to the 
ultimate tensile strength of the typical GRP ship hull 
material reported in Table 4. Figure 23 shows an equivalent 
loss of strength per decade for an initial ultimate tensile 
strength of 26,000 psi. The number of cycles to the extreme 
sea state slamming stress of 4700 psi which could be 
withstood by this material is calculated as 10?, a ya i ue 
unlikely to be. seen in any ship hulls. This treatment 

assumes that the material exhibits ideal mechanical 

properties and that any cutouts present have been adequately 
reinforced. However, this may not be the case since many 
unreinforced holes could exist in the ship hull. 

Stress concentration theory can be used to predict the 
failure point for a cutout that is exposed to cyclic 
loading. The predicted failure stress for a particular hole 
size is calculated by applying the stress concentration 
factor to the ultimate tensile strength and then reducing 
that value for each fatigue cycle. Any applied stress 
levels above these predicted stresses will cause failure. 
Figure 24 presents these data for two hole sizes as a 
function of stress and number of cycles. The plot shows 
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that -for slamming stresses o-f 4700 psi a 7.5 inch 
unrein-forced hole will -fail a-fter 10 cycles and that in seas 
without slamming 10,000 cycles will cause -failure. This 
model provides an alarming result that warrants a -further 
review o-f -fatigue using a -fracture mechanics approach. 

Linear elastic -fracture mechanics can be applied to 
■fatigue in notch sensitive woven -fabric and chopped strand 
mat composites by studying the propagation o-f macroscopic 
cracks [27]. A mode o-f ligamented crack growth, by which 
crack initiation is de-f i ned at a notch as the -first ligament 
■failure, is used in this analysis. This method is most 
easily de-fined by a -formula which predicts the number o-f 
cycles to crack initiation 



where 



log N = 




k i/Kq) 



(15) 



N = cycles to crack initiation 
ultimate tensile stress 
S = slope o-f S-N curve 

Kj= stress intensity -factor at the applied stress 
Kq= -fracture toughness 



This method can be used to predict the number o-f cycles to 
crack initiation -for the hull panels under study by using 
experimental values -for K Q and and a value o-f S -from the 

slope o-f Figure 23. These values are in close agreement 
with the values used in Reference 26 and are compared in 
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Table 11. The values of Kj -for each notch length and 
various applied stress levels are calculated by the method 
used in Appendix D. Figure 25 presents these data as a 
■function of applied stress and number of cycles to crack 
initiation for different notch lengths and shows that the 
number of cycles to failure decreases directly with 
increasing notch length. This plot is significant because 
once the first fibers are broken the crack propagation rate 
will increase significantly. This can be demonstrated by 
using the fatigue crack growth rate law [27] 

dc/cff^l = d/exp[ 2 . 3 <3 "f/S( 1 - Kj/k q >] (16) 

where d is the experimentally measured ligament width, which 
is assumed to be 0.125 in. based on measurements for a 
similar E-gl ass/pol yester material in Reference 26. Figure 
26 presents data from the crack growth rate law and shows 
the direct relationship between crack length and the crack 
growth rate. As the crack grows under a constant applied 
stress range, it accelerates rapidly. The analysis by 
fracture mechanics supports the findings of the stress 
concentration model which concluded that fatigue is a major 
determinant in the failure equation for a stressed GRP panel 
with cutouts and edge cracks. 

The use of high factors of safety in GRP ship design 
has precedent [28] and would be of benefit in reducing the 
chance of catastrophic failure in the hull. The findings in 
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this paper indicate that three conditions directly affect 
material strength and should be considered in determining 
factors of safety. Saltwater exposure, fatigue, and the 
presence of stress concentrations directly impact on the 
strength and are of primary concern. Figure 27 presents a 
simplified model to illustrate the components of a total 
factor of safety. The factors shown in Figure 27 include: 
(1) The effect of sea water on the strength, about 12. 9V. 
from Reference 26. <2) The effect of fatigue (in sea water) 

of 11. 67. per decade of cycles from Reference 26. The safety 
factor needed for various sea conditions, related to the 
expected number of cycles at a given maximum stress level 
can be assigned and used to account for the effects of 
fatigue. (3) The effect of stress concentration at 

unreinforced holes or other geometrical irregularities. The 
factors shown represent both the theoretical predicted 
stress concentration based on the measured ultimate tensile 
strength, and a lower factor which represents the actual 
experimental data. (4) Other factors which have not been 
addressed in this study, including statistical variations in 
material properties, flaws in fabrication, elevated 
temperatures, impact, etc.. 

The results shown in Figure 27 clearly indicate two 
areas which are critical in reducing the required safety 
factors. First, the reinforcement of holes and improved 
design and analysis procedures for other geometrical 
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irregularities. Second, the ability to operate in -Fatigue 
at higher stresses, either by improving the fatigue 
resistance or by raising the initial strength while 
maintaining a similar relative S-N curve slope. Graphite 

fiber reinforcement would improve the fatigue resistance (as 
well as the stiffness), while S-glass would raise the 
initial strength at a much lower cost than graphite. 



40 



CHAPTER 6 



CONCLUSIONS AND RECOMMENDATIONS 



The objective of this thesis was to help solve a 
practical, real world problem: to assess the effects o-f 

cutouts on GRP ship hull performance. The approach was to 
combine theoretical modeling and an experimental test 
program with specimens o-f appropriately large size. The 
scope precluded the analysis o-f many associated -factors, but 
the results resolve a specific portion o-f the technology gap 
in -fracture mechanics reported in Reference [13. 

Cone 1 us i ons 



The experiments conducted in this study provide data 
from which appropriate theoretical models for determining 
the failure strength of GRP hull panels can be selected. 
Test methods were developed for full scale tensile tests on 
16 large panels, 48 in. x 15 in., each with a central 
cutout. Ten of the panels had holes without edge cracks and 
six of the panels had holes with machined edge cracks. Two 
associated test methods, a standard tensile test and a 
conventional fracture toughness test, were conducted on a 
number of the remaining hull panels. All test data were 
then used for validation assessment of three theoretical 
mode 1 s . 
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The notch sensitivity was evaluated by compar i ng the 
ratio o-f the net strength o-f the remaining cross sectional 
area to the tensile strength -for panels with central holes, 
with and without machined edge notches. All ratios were 
much less than 1.0 indicating that the material is notch 
sensitive, so that a model which includes stress 
concentration e-f-fects is required. 

A stress concentration model was evaluated by comparing 
the experimental -failure stress in the 10 panels with simple 
holes to the predicted -failure stress, calculated using the 
measured ultimate tensile strength o-f the material and the 
orthotropic stress concentration -factor -for the particular 
hole size. This model appears to be valid based on the 
trend o-f the data, but an accurate prediction requires a 
local strength at the hole which is somewhat higher than 
that measured -from the standard tensile tests. 

For holes with edge notches, a linear elastic -fracture 
mechanics model was evaluated by comparing the 
experimentally applied stress at -failure with that predicted 
-from the -fracture toughness. The predicted values were 
close to the experimental results, but were slightly 
non-conservative. The -fracture mechanics model appears to 
be valid i -f notches o-f 1/4-inch or greater length are 
present at the edge o-f the hole. Discrepancies between 
measured and predicted -failure stresses may be due to the 
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isotropic analysis used and/or the large damage zone at the 
notch tip. 

The additional in-fluence of hull strength design 
practices, sea water environments, and fatigue cycling were 
evaluated by reference to other studies in the literature. 
Of particular importance was the effect of material fatigue 
when combined with the results for cutouts. To avoid the 
need for very high safety factors, cutouts in the size range 
studied here must be fully reinforced to eliminate their- 
stress concentrating effect. 

In summary, the need for high factors of safety and the 
requirement to reinforce large holes are two sound 
engineering practices which this research advises in using 
GRP for naval ship hulls. 

Recommendat i ons 

In concluding this thesis, the following recommendations 
are summarized to provide direction for continued research 
in the development of sound models for failure analysis in 
GRP hulls. 

- Development of a method by which the in-si tu strength 
at a hole boundary may be measured or predicted 

- Development of a series of guidelines for reinforcing 
hull penetrations in GRP ships, and assessment of the 
minimum hole size which requires reinforcement. 
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Investigation of the effect of different methods of 
machining cutouts on laminate strength. 



Full scale tensile tests on a large sandwich panel 
with cutouts 



Development of methods by which defects at cutout 
boundaries may be repaired. 



Development of a set of guidelines for inspecting hull 
penetrations. 

Development of a set of acoustic emission guideline 
data for sandwich panels with hull cutouts. Examination 
of the applicability of continuous in-service acoustic 
emission monitoring at critical stress points in an 
actual hull, based on guideline data. 

Consideration of a more fatigue resistant composite 
such as graphite/epoxy, or a higher initial strength 
material such as S-gl ass reinforcement, particularly 
in areas of high stresses. Further study is also 
needed of GRP in high cycle, low stress fatigue. 

Consideration of compressive loading effects. 
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TABLE 1 



Chief Advantages and Disadvantages of 
Glass Reinforced Plastics in Naval Applications [43 



Advantages 


Disadvantages 


High Strength/weight ratio 


- Low elastic modulus 


The ability to fabricate large, 
complex shapes in one piece 


- The differing elongations 
before failure of the two 



6RP constituents 

The ability to design directionally 

and therefore more economically - The highly anisotropic nature 



The comparative ease of forming 
complicated shapes 


of the material, in particular 
the relatively low 1LSS available 

- The combustibil ity of the material 


The resistance to corrosion 


- The relatively high material costs 


The resistance to marine 
biological attack 


- The difficulty or indeed the im- 


The ease of maintenance and 
repair 


practi cal i ty of checking the 
quality and ensuring the uniform- 
ity of large fabricated structures 


The fact that it is non-magnetic 


- The possible deterioration due to 
long term immersion in sea water, 


The variation of mechanical 


particularly under conditions of 



properties wi thin very wide limits high pressure and loading 



The fact that it does not warp, 
shrink or spl i t 


- The initial mold costs, for small 
numbers 


With increasing use, future costs 
are 1 ikely to ccme down 


- The fall off in appearance after a 
period in service 


The fact that it is radar and 
acoustically transparent 


-The inherent difficulty in joining 
parts together 


The fact that it has low thermal 
conducti vi ty 




The high degree of surface finish 
needs to be produced only on the 
mold surface 





The ability to use sandwich material 
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TABLE 2 



Forecast o-f Demand -for Rein-forced Plastics. 1975-1985 [53 



Mill ions of lb 


1985 


1984 


1983 


1982 


1981 


19®) 


1975 


Transportation, land 


570 


540 


458 


359 


445 


414 


245 


Construction 


443 


430 


400 


312 


309 


287 


175 


Anticorrosion 


324 


310 


288 


235 


275 


252 


143 


Marine 


325 


309 


274 


230 


290 


275 


285 


Electrical 


205 


189 


170 


140 


178 


142 


82 


Consumer goods 


149 


143 


128 


84 


110 


103 


44 


Appliance 


130 


123 


104 


82 


112 


104 


44 


Aircraft 


35 


29 


25 


22 


28 


25 


24 


Other 


84 


80 


72 


44 


74 


70 


53 


TOTAL 


2249 


2153 


1923 


1528 


1821 


1494 


1175 



Annual Chanoe % 


1984-1985 


1983-1984 


1980-1985 


1975-1985 


Transportation, land 


5.4 


17.9 


4.5 


8.0 


Construction 


3.0 


7.5 


9.1 


9.7 


Anticorrosion 


5.2 


7.4 


5.3 


7.2 


Marine 


5.2 


12.0 


3.4 


1.3 


Electrical 


8.5 


11.2 


4.8 


9.4 


Consumer goods 


4.2 


11.7 


7.7 


8.8 


Appl iance 


5.7 


14.0 


4.4 


7.3 


Aircraft 


20.7 


14.0 


7.0 


3.8 


Other 


7.5 


11.1 


4.2 


5.0 


TOTAL 


5.4% 


12.07. 


4.0% 


4.8% 
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TABLE 3 



Existino Mine War-fare Ships with GRP Hulls [10] 



Country 


Class 


Length 

(ft) 


Displ . 
(tons) 


Entered 

Service 


Number 

Planned 


Notes 


Austral ia 


Catamaran 


101 


160 


Mid-80's 


6 


1 


Belgium 


Tripartite 


154 


544 


80-87 


10 


2 


France 


Tripartite 


154 


544 


80-87 


15 


2 


Italy 


Lerici 


163 


500 


1983 


10 


3 


Malayasia 


Hah am i ru 


163 


500 


1984 


4 


3 


Netherlands 


Tripartite 


154 


544 


80-87 


15 


2 


Sweden 


MSO 


155 


340 


1984 


6 


1 


Britain 


Wilton 


152 


450 


1973 


1 


2 


Britain 


Hunt 


197 


725 


1980 


12 


2 


USSR 


Zhenya 


140 


300 


1970 


3 


- 


USSR 


Andryusha 


147 


360 


1970 


2 


5 


USSR 


Sonya 


160 


460 


1973 


42 


4 


USA 


MSH 


150 


434 


1987 


3 


1 



Notes: 

1. GRP foam core sandwich 

2. GRP single skin, stiffened 

3. GRP single skin, monocoque 

4. Wood with GRP sheath 

5. GRP structure, unconfirmed 
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TABLE 4 



GRP Naval Ship Hull Material Properties [22] 



Hull materials combination of marine grade woven 

roving and chopped strand mat in a 
marine grade isophthalic polyester 
resin 



Glass content: 45X + 3Z. by weight 



Woven Roving: 18oz./yd2 

Chopped Strand: Poz./yd^ 

Minimum Mechanical Properties: 

24,200 psi 
1 .7 x 10* psi 

38.000 psi 
1 .7 x 10* psi 
25 , OOOps i 
1 .9 x 10* psi 

14.000 psi 
0.45 x 10* psi 

1.0 x 10^ psi 

0 . 1 ? 



tensile strength 
tens i 1 e modu 1 us 
■flexural strength 
■flexural modulus 
compressive strength 
compressive modulus 
shear strength 
shear modulus 
I LSS 

Poisson's ratio 
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TABLE 5 



Maximum Orthotropic Stress Concentration 
Factors and Predicted Applied Failure 
Stress as a Function o-f Hole Diameter 



Hole Diameter(in) 



K 



Tor thotrop i c 



■fai 1 ure < ps i ) 



3.75 


3.727 


6975 


5.00 


3.998 


6503 


7.00 


4.732 


5494 


7.50 


4.954 


5248 




*assum i ng "^ts = 26000 psi 
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TABLE 6 



Number o-f 
Samp 1 es 



Full Size GRP Test Panel 
Charac ter i st i cs 



Hoi e 

D i ame ter 


Crack 

Length 


Mach i n i ng 
Central Cutout 


7.50 





1 


5.00 


- 


1 


5.00 


- 


2 


3.75 


— 


1 


5.00 


1 .00 


1 


5.00 


0.50 


1 


5.00 


0.25 


1 


Carbide tipped parting 
in a -flycutter 
Arbor mounted hole saw 


tool mounted 
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TABLE 7 



Characteristics o-f and Tensile Tests 
Results -for Large GRP Panels 



Panel 

Number 



Hole Crack 

Diameter(in) Length(in) 



Fai 1 ure 
Load(lb-f) 



1 


5.00 


- 


44560 


14259 


2 


5.00 


- 


46990 


15037 


3 


5.00 


- 


47540 


15213 


4 


5.00 


- 


47260 


15123 


20 


5.00 


- 


48240 


15437 


21 


5.00 


- 


50530 


16170 


12 


3.75 


- 


52470 


14925 


13 


3.75 


- 


52870 


15039 


5 


7.00 


- 


41260 


16504 


14 


7.50 


- 


36940 


15761 


15 


7.50 


— 


36300 


15488 


6 


5.00 


1 .0 


27350 


10940 


7 


5.00 


1 .0 


25080 


10032 


8 


5.00 


0.5 


29980 


10659 


9 


5.00 


0.5 


29320 


10425 


1 0 


5.00 


0.25 


35030 


11799 


1 1 


5.00 


0.25 


30980 


10435 


1? 


Ful 1 




23620 


27485 


18 


Size 




22360 


26937 


19 


Tens i 1 e 




18860 


23552 
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TABLE 8 



Conventional Tensile Test Results 



Samp 1 e 


Spec i men 


Failure 


Number 


Orientation 


StressCpsi > 


22 


Vertical 


27554 


23 


Vertical 


24960 


24 


Vertical 


30933 


25 


Hor i zontal 


32000 


26 


Hor i zon tal 


35520 


27 


Hor i zon tal 


35413 



Young 7 s 
Modul us<psi ) 



1 .62 x 106 
1.84 x 10* 

2.47 x 10* 
2.42 x 10* 
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TABLE 9 



Fracture Touohness Test Results 



Sample Specimen Propagation Crack Isotropic Orthotropic 
Number Orient. Load <lbf> Length Kq Kq 



28 


Ver t i cal 


1725 


1 .484 


26.2 


22.9 






1525 


2.359 


31 .2 


29.2 






1 175 


2.625 


25.9 


24.8 


29 


Vertical 


1850 


1.812 


31 .8 


28.4 






1550 


2.219 


30.4 


27.7 






1375 


2.469 


29.0 


27.9 


30 


Vertical 


1875 


1 .562 


29.3 


26.4 






1675 


1 .844 


29.0 


26.2 






1325 


2.250 


26.2 


24.1 


31 


Hor i zon tal 


2175 


1 .360 


31 .4 


27.8 


32 


Hor i zon tal 


2275 


1 .312 


32.2 


28.4 


33 


Hor i zon tal 


2175 


1 .375 


31 .6 


27.8 
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TABLE 10 



Stress Concentration Model Predicted 
Stress Versus Measured Failure Stress 



Hoi e 

D i ame ter < i n ) 



pr ed 



..fa i 1 



3.75 


6976 


11236 


5.00 


6 514 


10137 


7.00 


5489 


8802 


7.50 


5246 


7812 



Ratio 



0.620 

0.642 

0.623 

0.671 
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TABLE 1 1 



Dynamic Fatigue Properties -for 
Woven Rovino/Mat Laminates in Salt Ulater 



Re-ference 1261 Table 4 



Measured (|<si) 


31 .0 


26.0 


Extrapolated <ksi) 


27.0 


22.6 


S (ksi/decade) 


3.14 


3.02 


Kq (ksi in ) 

Ligament Width, d (in.) 


24.8 


25.9 


.125 


.125 
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FIGURE 1 



Geometry -for Notch Insensitivity Model 



P 




P 



5 ? 



FIGURE 2 



P 



Isotropic Finite Width Plate 




P 



P = uni-form tensile stress 
a = hole radius 
b = plate hal-f width 

X = a/b 
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FIGURE 3 






Orthotropic Infinite Plate 



ii 




n 



a 




v 



t 







t 
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FIGURE A 



Finite Plate 
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FIGURE 5 



Stress Intensity Factor Graph for a 
Rectangular Plate with Two Edge 
Cracks at a Circular Hole (From [171) 
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